DNS of premixed turbulent V-flame: coupling spectral and finite difference methods

To allow for a reliable examination of the interaction between velocity fluctuations, acoustics and combustion, a novel numerical procedure is discussed in which a spectral solution of the Navier–Stokes equations is directly associated to a high-order finite difference fully compressible DNS solver (sixth order PADE). Using this combination of high-order solvers with accurate boundary conditions, simulations have been performed where a turbulent premixed V-shape flame develops in grid turbulence. In the light of the DNS results, a sub-model for premixed turbulent combustion is analyzed. To cite this article: R. Hauguel et al., C. R. Mecanique 333 (2005).     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

爆炸地震波评价地面毁伤效应的实验研究

通过实验探索爆炸地震波评价爆源对地面冲击波毁伤效应的可能性 ,结果表明 ,在相同的地层状态条件下和一定的距离范围内 ,用爆炸地震波评价爆源对地面的毁伤效应是可行的。这在很大程度上弥补了利用空气冲击波压力测试和模拟目标毁伤评价燃料空气炸药 (FAE)爆源爆炸空气冲击波威力的不足。     

冲击伤复合高速破片致伤的损伤特点和冲击伤防护研究

探讨了冲击伤复合高速破片致伤的损伤特点及复合材料对冲击伤的防护效果。结果表明，肢体高速破片致伤对肺冲击伤具有加重作用，加重效应主要发生在原发冲击伤部位，加重效应的机理可能与压力波引起的血流扰动有关。利用复合材料防护可减轻肺冲击伤程度，降低动物的死亡率。     

Particle behavior in homogeneous isotropic turbulence

Direct numerical simulations were conducted to investigate the behavior of heavy particles in homogeneous isotropic turbulence. The present study focused on the effect of particle inertia and drift on the autocorrelations of the particle velocity and the fluid seen by particles and the dispersion characteristics of particles. The Lagrangian integral time scale of particles monotonically increased as the magnitude of the particle response time increased, while that of the fluid seen by particles remained relatively constant; it reached a maximum when the particle response time was close to the Kolmolgorov time scale of the flow. Particle dispersion increased as the particle inertia increased for small particles, while for larger particles, it decreased as particle inertia increased; particle eddy diffusion coefficient was maximal, and greater than that of the fluid by about 30%, at the preferential concentration. The concentration field of the particles with τp/τk≈1.0 showed that particles tend to collect in regions of low vorticity (high strain) due to preferential concentration. As the drift velocity of a particle is increased it crosses the paths of fluid elements more rapidly and will tend to lose correlation with its previous velocity faster than a fluid element will. And the correlation of particle velocities along the drift direction is more persistent than that perpendicular to the direction of drift. Simulations also showed that the continuity effect and the crossing-trajectory effect are weakened for particles with infinite inertia.     

Direct modeling of flow of FENE fluids

Direct simulations of macromolecular fluids are carried out for flows between parallel plates and in expanding and contracting channels. The macromolecules are modeled as FENE dumbbells with soft disks or Lennard-Jones dumbbell-dumbbell interactions. The results are presented in terms of profiles and contour plots of velocity, pressure, temperature, density, and flow fields. In addition the data for potential energy, shear stress, and the normal components of the stress tensor are collected. In general, an excellent agreement is found between the simulated profiles and the well-known flow structures, such as flow separation and formation of viscous eddies, indicating that micro-hydrodynamics is a viable tool in linking macroscopic phenomena with the underlying physical mechanisms. The simulations are performed in the Newtonian regime, for medium-size systems comprising up to 3888 dumbbells. This number is sufficiently large to control boundary and particle number effects. The flow is induced by gravity. The traditional stochastic (thermal) and periodic boundary conditions are employed. Also, diffusive boundary conditions, which could include a stagnant fluid layer and repulsive potential walls, are developed. The scaling problems, which are related to the application of a large external force in a microscopic system (of the size of the order 100 Å), result in extreme pressure and temperature gradients. In addition, the viscosity and thermal conductivity coefficients obtained from velocity and temperature profiles of the channel flow are presented. These results are confirmed independently from modeling of Couette flow by the SLLOD equations of motion and from the Evans algorithm for thermal conductivity.     

爆炸荷载作用下钢筋混凝土梁的动力响应及破坏形态分析

对基于Timoshenko梁理论建立的非线性动力有限元法作了改进。根据压区理论得到混凝土的平均剪应力和平均剪应变关系,建立了能反映箍筋的抗剪作用的材料模型;此外,对结构在爆炸荷载作用下可能出现的各种响应现象进行了描述,以准确地预测梁破坏时不同位置截面上钢筋和混凝土的受力、变形及破坏情况。应用改进的材料模型,对爆炸荷载作用下的五个钢筋混凝土试验梁的动力响应和破坏形态进行了数值模拟,结果表明,该数值方法能较好地模拟钢筋混凝土梁的弯曲、弯剪和剪切等破坏形态。     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.